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1. Introduction 

The electromagnetic form factors of the proton and neutron arc quantities of continuous inter¬ 
est both experimentally and theoretically. Electromagnetic scattering of nucleons reveals properties 
such as their charge and magnetism, with insight on the distribution of charge within the nucleon. 
More recently, precision experiments have revealed a surprising discrepancy in the charge radius 
of the proton. Namely, the proton radius when measured recently via the Lamb shift of muonic hy¬ 
drogen [TJ has a value that is smaller by five standard deviations as compared to experiments using 
electron scattering and hydrogen Lamb shift |2]. Lattice QCD can provide insight on the origin of 
this discrepancy by evaluating from first principles the QCD contribution to the proton charge. So 
far, although consistent across discretisation schemes, lattice simulations at heavier than physical 
pion masses have underestimated the proton charge With simulations directly at the physical 
pion mass now becoming available, chiral extrapolations are no longer required, thus eliminating 
one source of systematic uncertainty and allowing direct contact with experiment. In this proceed¬ 
ings contribution we present preliminary results of the nucleon electromagnetic form factors using 
three ensembles of twisted mass fermion (TML) configurations; two Nf = 2 + 1 + 1 ensembles at 
pion mass 373 MeV and 210 MeV and one Nf = 2 ensemble at the near physical pion mass value 
of 130 MeV. 

2. Lattice Setup and Methods 

The electromagnetic matrix element of the nucleon can be decomposed into the Dirac (FQ and 
Pauli (Ff) form factors: 

(N(p'J)\f\N(p, S )) = ^ = 7^,V) + ^F2(T>, 

( 2 . 1 ) 

with m-N the nucleon mass, p (s) and // (s') the initial and final momentum (spin) of the nucleon, 
u are fermion spinors and q = p' — p is the momentum transfer. Alternatively one can define 
the so-called Sachs electric and magnetic form factors: GE(q 2 ) = F\(q 2 ) + Ffiq 2 ) and Gm(t) = 
F\(q 2 ) + TFjiq 2 ) with T = jsp— 1 - The slope of the form factors at zero momentum defines the 
relevant radii, namely the Dirac and Pauli radii via: (rf) = — j| j^\q 2 =o and similarly for the electric 
and magnetic radii (r|) and (r]^). 

On the lattice, one calculates an appropriate three-point correlation function: 

G^r-q-t s , ti ) = (2.2) 

XjX,- 

with Xn the nucleon interpolating operators, x s , x; and xo the final (sink), insertion and initial 
(source) coordinates and the local electromagnetic current: /' (x) = hl(x)Y‘‘u(x) — \ ] d(x)y- > d(x). 
In this work we use the lattice conserved electromagnetic current (3), and thus lattice results need 
no renormalisation. In our setup we fix p' = 0 and cany out the sum over x v using a sequential 
inversion through the sink. Lor this we require t s — to and the projection matrix T to be set before 
the sequential inversion. The choices F 4 = |(1 + 74 ) and T^. = Y.j^'^YsYj isolate Gf and Gm 
respectively. 
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Figure 1: Diagrams contributing to the nucleon three-point correlation function, connected (left) 
and disconnected (right). 


The three-point correlation function has quark-connected and quark-disconnected contribu¬ 
tions as shown in Fig. [T] Here we compute the isovector (V) and isoscalar (,S) combinations of 
the operator: j ,Ll s (x) = u{x)y p u(x) + cl(x)Y‘‘ d(x). Assuming SU(2) flavour symmetry the discon- 
nected contribution cancels for the isovector current, but does not in the isoscalar combination. 
Such disconnected contributions are exceptionally difficult to compute on the lattice due to the 
increased susceptibility they exhibit to statistical fluctuations. However, recent results using new 
techniques have shown that at low momentum transfer they are small, of order of 1 %, compared 
to the connected, at pion masses of around 300 to 400 MeV (6}[7j. From the isoscalar and isovector 
combinations we can obtain the proton and neutron form factors through the linear combinations: 
F p -F n = F u -F d and F p + F n = \(F U + F d ). 
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Figure 2: The lattice spacing is set using the nucleon mass. 
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Vf = 2 + 1 +1, csw 

= 0 

1.95 

32 3 x64 

0.082 

0.373 

2.10 

48 3 x96 

0.064 

0.210 

Nf = 2, csw 

= 1.57551 

2.10 

48 3 x96 

0.094 

0.130 


Table 1: The twisted mass fermion 
ensembles used in this work. 


In this work we present three ensembles of TMF configurations; two N\ = 2+1 + 1 ensembles 
with m K = 373 MeV (referred to as “B55”) and 210 MeV (referred to as “D15”) and one Nf = 2 
ensemble with m n = 130 MeV (8), referred to as the physical point ensemble. More details arc 
given in Table [T] The lattice spacings are determined by fitting the nucleon masses of 18 TMF 

icavy baryon chiral perturbation theory expansion: tun = m° N — 
We require the curve to reproduce the physical nucleon mass 


ensembles using the leading order 1 

r\ ^ Q 

4 ci»7 ^. — * A f2 m K , as shown in Fig. 


16 71 fl' 


and allow the spacings to vary as fit parameters. This yields the spacings listed in Table [T] for the 
ensembles used here. 


For extracting the matrix element from the lattice data, we form an appropriate ratio of three- 
point to two-point functions, which cancels unknown overlaps and exponential energy factors such 
that in the limit of large time separations: — to^$> and t s — tj ^$>, the contributions from excited 
nucleon states are damped out, leaving a time-independent ratio (referred to as the plateau region)-. 
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Figure 3: Extraction of the required matrix element is done via fitting to the plateau method (left) 
or the summation method (right). In the left panel we compare the result of the summation method 
(black line) with that of the plateau (yellow line) for the B55 ensemble. 


R v M (ti,t s \ k 


2 \ 


f.—«o» 


> G^(k 2 )[1 + o(e~ AM ^ ts ~ ti ^ ,g-A^OOfr— f o)] 


(2.3) 


as illustrated for the isovector magnetic form factor in Fig. [3] AM (A E) is the mass (energy) 
gap of the nucleon ground to nucleon first exited state. With fixed sink time-slice, we vary the 
insertion time-slice to identify and tit to the plateau region. This will be referred to as the plateau 
method. To ensure ground state dominance, we invert for multiple sink-source separations t s — to, 
with the incremental EigCG algorithm as an efficient multiple right-hand-side solver Q. With the 
availability of multiple sink-source separations, we can also sum over the insertion time-slice t\ to 
obtained the summed ratio as a function of t s : 


£<(t i) 4;k 2 )-^>C + G^(k 2 )4[l + O(e- AM ^- f0 ), e - AE «fe- f °))]. (2.4) 

u 

This allows an alternative extraction of the form factor via a two parameter linear fit, the so-called 
summation method. Compared to the plateau method, the advantage is that excited states decay 
with a larger exponential suppression factor determined by t s — to rather than t s — U or f, — to, while 
the disadvantage is that a two parameter fit is required. In Fig. [3] we show an example of fits to the 
summed ratio for momenta up to k 2 = 6. 


3. Results 

In Fig. [4] we show results for the B55 ensemble for which we have carried out inversions 
at seven sink-source separations for 1,200 configurations. We show the electric and magnetic 
isovector Sachs form factors obtained via the plateau method and compare them to the summation 
method to identify excited state effects. As the separation increases from 0.5 to 1.4 fm the form 
factors become steeper, yielding larger electric and magnetic radii. However these results are still 
in tension with experiment, shown by the curve which represents J. Kelly’s parameterisation of the 
experimental results. 
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Figure 4: The electric (left) and magnetic (right) nucleon form factors versus the momentum trans¬ 
fer squared, for B55 for various values of t s (coloured symbols) and for the summation method (sm, 
black asterisk). The solid line is J. Kelly’s parameterisation to the experimental data. 


The radii are determined from the derivative of the form factors at Q 2 = 0. Fitting the form 
factors to a dipole form Fj(Q 2 ) = F,(0)/[1 + Q 2 /M ~] 2 the radii are determined from the dipole 
mass via (rj) = 12 /Mf with i = 1,2, and similarly for the Sachs form factors to obtain the electric 
and magnetic radii. For the Dirac and electric form factors F\ (0) = Ge( 0) = 1 since we are using 
the conserved current and thus only the mass is a fit parameter. Fj(0 ) is fitted for and yields the 
anomalous magnetic moment of the nucleon. 
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Figure 5: The Dirac and Pauli isovector radii (left and centre), and the isovector magnetic moment 
(right) for various t s and for the summation method. The notation is the same as in Fig. |4| 


In Fig. [5] we show the isovector Dirac and Pauli radii and the isovector anomalous magnetic 
moment for the B55 ensemble. Results extracted via fits to the dipole form as a function of the 
sink-source separations, as well as the value extracted from the summation method. These results 
corroborate the conclusion that as t s increases the radii become larger. We also observe that con¬ 
tamination from excited states are suppressed after a sink-source separation of around 1.2-1.3 fm, 
indicated by the agreement between plateau and summation methods beyond this distance. 

Our results using all three ensembles, including the physical point ensemble, are shown in 
Fig. [6j for the isovector electric and magnetic Sachs form factors. A single sink-source separation 
of 1.2 fm is available for the m K =2 10 MeV ensemble, while for the physical point we have sink- 
source separations at 0.94, 1.13 and 1.32 fm, the largest of which is shown in Fig. [6| for around 
1,400 configurations. We observe a steeper form factor as the pion mass is reduced, although more 
statistics are required at the physical point for a definite conclusion. In the same figure the electric 
and magnetic radii are also shown as a function of the pion mass, confirming the tendency towards 
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Q 2 [GeV 2 ] Q 2 [GeV 2 ] [GeV] 


Figure 6: The isovector electric (upper left) and magnetic (lower left) form factors for the three 
TMF ensembles of this. The associated radii are shown in the right panel as a function of the pion 
mass. 


larger values as the pion mass is decreased. 

We compare our results at the physical point with those of a recent calculation using clover 
fermions at similar pion mass in Fig. |T| fTOj . We see a consistency in the extracted form factors 
between the two formulations. The isoscalar Pauli form factor F^, not shown in Fig. [ 7 ] is found 
consistent with zero in both cases. 



Q [GeV 2 ) Q [GeV 2 ) Q [GeV 2 ) 

Figure 7: The isovector Dirac (left) and Pauli (centre), and the isoscalar Dirac (right) form factors 


of this work, compared to results from the LHPC at m n = 149 MeV [ 101. 


In Fig. [8] we show the isovector Dirac and Pauli radii as a function of the pion mass, for the 
TMF ensembles analysed in this work and the clover ensembles analysed in Ref. 1101. We also 
show the PDG values of these quantities and, for the case of (r \) v , we include the recent result 
obtained from muonium Lamb shift measurements. As the pion mass decreases the effect of the 
excited states seem more pronounced indicated by the spread of the open symbols. However, the 
errors also become larger and a study with increased accuracy is required. Both clover and TMF 
results are in agreement, and seem to be converging towards the experimental values, although 
errors are large and need to be drastically reduced to make contact with experiment. 


4. Conclusions 

We have presented preliminary results for the electromagnetic form factors of the nucleon 
including a lattice QCD ensemble with a pion mass set to its physical value. A lattice at a heavier 
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m„ [GeV] m n [GeV] 

Figure 8: The isovector Dirac (left) and Pauli (right) radii as a function of the pion mass. The 
diamonds are from this work while the squares are from the LHPC. Filled symbols are results 
using the summation method, while open symbols are from the plateau method. 

pion mass of 373 MeV was used for a thorough study of excited state effects, indicating that 
a three-point function sink-source separation of ~1.3 fm is sufficient for ensuring ground state 
dominance. Our results at the physical point are in agreement with recent results using clover 
fermions at a similar pion mass. Within large statistical errors, our results are either in agreement 
or tend towards the experimental values as the pion mass is reduced. In order to make contact with 
experiment, a 2% error in the radii is required, which amounts to more than a 10-fold increase in 
statistics. Deflation algorithms in combination with all-mode-averaging [111] are currently under 
investigation to efficiently achieve this. 
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